Self-similar chain conformations in polymer gels 
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We use molecular dynamics simulations to study the swelling of randomly end-cross-linked polymer net- 
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works in good solvent conditions. We find that the equilibrium degree of swelling saturates at Q eq « N e 
for mean strand lengths N s exceeding the melt entanglement length N e . The internal structure of the network 
strands in the swollen state is characterized by a new exponent v — 0.72±0.02. Our findings are in contradiction 
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to de Gennes' c -theorem, which predicts Q cq ~ N s and v = 0.588. We present a simple Flory argument 
for a self-similar structure of mutually interpenetrating network strands, which yields v = 7/10 and otherwise 



recovers the classical Flory-Rehner theory. In particular, Q eq 
PACS Numbers: 61.41+e,82.70.Gg,64.75.+g 
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if N e is used as effective strand length. 



Polymer gels [jl|-^] are soft solids governed by a complex 
interplay of the elasticity of the polymer network and the poly- 
mer/solvent interaction. They are sensitive to the preparation 
conditions and can undergo large volume changes in response 
to small variations of a control parameter such as tempera- 
ture, solvent composition, pH or salt concentration. In this 
letter we reexamine a classical but still controversial problem 
of polymer physics |jl],^|], the equilibrium swelling of a piece 
of rubber in good solvent. 

Experimentally gels have been studied extensively by com- 
bining thermodynamic and rheological investigations with 
neutron or light scattering . Here we use computer sim- 
ulations [ffl-|l3[l, since they offer some advantages in the ac- 
cess to and the control over microscopic details of the network 
structure. We concentrate on the role of entanglements in lim- 
iting the swelling process and, in particular, the structure of 
the network strands in the swollen gel. Questions relating to 
the structural heterogeneity in our gels and the butterfly ef- 
fect [p|Jl4l] will be addressed in a future publication. 

There are several theories for the swelling of polymer net- 
works prepared from a melt of linear precursor chains. In the 
dry state of preparation, the network strands have Gaussian 
statistics, i.e. the mean square end-to-end distance is related to 
the average length, N s , by {r 2 ) dry oc b 2 N s 2u , where v = 1/2 
and b is the monomer radius. The same relation also holds 
for all internal distances, leading to the characteristic struc- 
ture factor S(q) ~ q~ x / u for the scattering at wave vector q 
from a fractal object. 

The classical Flory-Rehner theory [|l]] writes the gel free en- 
ergy F as a sum of two independent terms: a free energy of 
mixing with the solvent (favoring swelling and estimated from 
the Flory-Huggins theory of semi-dilute solutions of linear 
polymers) and an elastic free energy (due to the affine stretch- 
ing of the network strands which are treated as Gaussian, 
concentration-independent, linear entropic springs). Mini- 
mizing F yields Q eq oc N s 3 ^° for the equilibrium degree 
of swelling. The Flory-Rehner theory implies that the struc- 
ture factor of long paths through the network is of the form 
S(q) ~ q~ 2 both locally, where the chains are unperturbed, 
and on large scales, where they deform affinely ((r 2 ) eq oc 



{r 2 ) dry Q% 3 ) with the outer dimensions of the sample. The 
stretching should be visible in the crossover region around 
q ss 2ir/{bN s 1/2 ) with S(q) ~ q~ l . 

More recent treatments are based on the scaling theory of 
semi-dilute solutions of linear polymers [Hj . Locally, inside of 
so called "blobs", the chains behave as isolated, self-avoiding 
walks with v w 3/5. A blob containing g monomers has a 
typical diameter £ oc bg v — bQ 3 ^ 4 . Chains with N 3> g 
can again be regarded as ideal, however, with a renormalized 
chain length N/g and a renormalized monomer size £. Quite 
interestingly, refining the Flory-Rehner ansatz along these 



lines [ir&|19|J14|] recovers the classical result Q e 
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Such models imply for the structure factor a crossover from 
S(q) ~ q~ 5/3 to S(q) ~ q^ 1 atgw 2ir/£ eq , where the blob 
diameter at equilibrium swelling is given by £ eg oc bQ 3 q 4 and 
much smaller than the strand extension. 

An open point is the length scale on which the systems be- 
gin to deform affinely with the macroscopic strain. The two 
theories mentioned above consider linear entropic springs, 
where due to the global connectivity (disregarding fluctuation 
effects) this length scale is given by the strand size. A drasti- 
cally different point of view has been advanced by de Gennes, 
who argues that the swelling is limited by the local connectiv- 
ity, which only begins to be felt at the overlap concentration 
c* oc N s /(bN s 1 ') 3 of a semi-dilute solution of linear poly- 
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mers of average length N s , corresponding to Q e 
As a motivation for his c*-theorem [Q], de Gennes consid- 
ers crosslinking in dilute solution, but postulates that the same 
results also hold for swelling of networks prepared by cross- 
linking dense melts. The c*-theorem predicts S(q) ~ g~ 5 / 3 
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for q > 2-K/(bN s ) as well as unusual elastic properties due 
to the non-linear elasticity of the network strands Jl2[|T5| ]. 

Both, the Flory-Rehner theory |jl|] and the c* -theorem [0] 
are supported by part of the experimental evidence [^, |T^|T7| ]. 
While the results are very sensitive to the details of the prepa- 
ration process, it seems well confirmed [|l7[| that highly cross- 
linked networks behave according to Flory's prediction Q eq oc 
N 3/5 . On the other hand, SANS experiments [|,|7|] and com- 
puter simulations [[l(]] of lightly cross-linked gels show the 
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weak dependence of the strand extensions on the degree of 
swelling predicted by the c* -theorem. To our knowledge, all 
SANS studies have concentrated on the low q Guinier regime 
(i.e. the radius of gyration of the strands) and no particular 
attention has been paid to the local chain structure. 

As in earlier investigations of polymer melts and networks 
|[2Cj-|24|], we used a coarse-grained polymer model where 
beads interacting via a truncated, purely repulsive Lennard- 
Jones (LJ) potential are connected by FENE springs. With e, 
cr and r as the LJ units of energy, length and time, the equa- 
tions of motion were integrated by a velocity-Verlet algorithm 
with a weak local coupling to a heat bath at ksT — le. The 
potentials were parametrized in such a way that chains were 
effectively uncrossable, i.e. the network topology was con- 
served for all times. In our studies we did not simulate the 
solvent explicitly, but rather used vacuum which can be con- 
sidered as a perfect solvent for our purely repulsive (athermal) 
network chains. The relevant length and time scales for chains 
in a melt are the average bond length, \J (I 2 ) = 0.965(5)<r, the 
mean-square end-to-end distance {r 2 )(N)dry = 1.74(2)l 2 N 
[{2"o|], the melt entanglement length, N e = 33(2) monomers, 
and the Rouse time TR OUSe (N) ~ 1.35rA 2 p5|]. In dilute so- 
lutions, single chains adopt self-avoiding conformations with 
(r 2 )(N) « 1.8/ 2 iV 3 / 5 . 

Using this model, it is possible to study different network 
structures including randomly cross-linked, randomly end- 
cross-linked [21 22j] and end-linked melts [|23|] as well as net- 
works with the regular connectivity of a crystal lattice p4|]. 
Here we investigate end-cross-linked model networks created 
from an equilibrated monodisperse melt with M precursor 
chains of length A at a melt-like density pd ry = 0.85cr~ 3 
by connecting the end monomers of the chains to a ran- 
domly chosen adjacent monomer of a different chain. This 
method yields defect-free tri-functional systems with an ex- 
ponential distribution of strand lengths N s with an average of 
N s = N/3. The Gaussian statistics of the strands remains 
unperturbed after crosslinking [ [?lp6| ]. The systems studied 
range from M/N = 3200/25 (i.e. the average strand size 
N s = 8.3) up to M/N = 500/700 (N s = 233), some sys- 
tems being as large as M ■ N = 5 • 10 5 . All simulations used 
periodic boundary conditions in a cubic box and were per- 
formed at constant volume. Starting from Vd ry = MN/ pdry, 
the size of the simulation box was increased in small steps 
alternating with equilibration periods of at least 5 entangle- 
ment times Tn{N e ) w 1400r. The isotropic pressure P was 
obtained from the microscopic virial tensor and the condi- 
tion P eq = was used to define equilibrium swelling with 
Qeq = V e q/Vdry Tests with a part of the networks using 
open boundaries did not show any significant changes of the 
results. 

We investigated the equilibrium swelling of our model net- 
works as a function of the average strand length N s . Fig. [j] 
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shows QeqNe^ 5 as a function of the average strand length 
(N e /N s )~ 3 / 5 . Our results for short strands are compatible 
with the Flory-Rehner [[!]] prediction Q 



A„ 3/5 , but do 
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( N e /N s! 

FIG. 1. Strand length N a dependence of the equilibrium de- 
gree of swelling Q eq - The straight line going through zero represents 
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Flory's prediction Q eq oc N s . The melt entanglement length N e 
was used to normalize the axis in order to show that deviations from 
Flory's theory occur around N s ~ N e and that the asymptotic value 
Qeq{N s — > 00) is of the order of N^ 5 . 



not allow for an independent determination of the exponent. 
In contradiction to this theory we observe a saturation of the 
equilibrium swelling degree for large N s . The crossover oc- 
curs for N s « N e . The extrapolated maximal degree of 
equilibrium swelling Q ma ,x(N s — > 00) = 6.8(3) is close to 

the swelling degree of an ideal Flory-gel with average strand 
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length N e : 1.15 ■ N e — 9.5, where the prefactoris empiri- 
cally obtained from the slope of the straight line in Fig. [j]. In 
contrast, the corresponding estimate based on the c* -theorem, 
Q eq ~ 6 3 /er 3 A e 4/5 » 36, is clearly too high {b = 1.3a is the 
stastical segment length in good solution). Our interpretation 
is that to a first approximation entanglements act as chemi- 
cal crosslinks in limiting the swelling of polymer networks. 
The situation is analoguos to an "olympic gel" [^] of topolog- 
ically linked ring polymers. In contrast to solutions of linear 
polymers, systems containing trapped entanglements cannot 
be arbitrarily diluted. 

The chain conformations at equilibrium swelling are best 
characterized by their structure factor S(q). Fig. |5| shows 
S(q) of the precursor chains within the network for our most 
weakly crosslinked N — 700 sample. We have chosen the 
Kratky-representation {q 2 S(q) vs. q) to show the deviation 
from the Gaussian case (S(q) oc q~ 2 ) more clearly. The 
observed power law form S(q) oc q~ x l v is characteristic of 
fractals and common in polymeric systems. However, the ob- 
served exponent v = 0.72(2) is unexpected. Furthermore, 
the fractal structure is observed for a g-range of 2er < =^ < 
15.5fi ~ bN®- 72 , suggesting that the mean extension of the 
effective strands of length N e is the only relevant length scale 
in the problem. For smaller q we see the onset of the expected 
scattering of a Gaussian chain consisting of randomly oriented 
parts of length bN®- 72 . Our precursor chains (even N = 700) 
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FIG. 2. Single chain structure function in networks at equilib- 
rium swelling in Kratky representation. The straight line corresponds 
to a power law q 2 ~ (u = 0.72). The figure contains scattering 
data for the precursor chains (•) of length N = 700 (N s — 233) and 
for network strands of lengths N s = 10 (O), N s = 40 (A), N s = 70 
(0), N s = 100 (□) within a different network with precursor chains 
of length N= 100 (N s =33). 
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FIG. 3. Log-log plot of the mean-square end-to-end distance 
{r 2 } eq of the individual network strands within a single network 
(N a — 33) at equilibrium swelling Q — 5.8 versus strand length 
N s . The straight line corresponds to a power law (r 2 ^ ~ N%" 
with v — 0.72. The data are normalized to an affine deformation 



are to short to see it clearly developed. 

Since the scattering from the precursor chains could be 
affected by polydispersity effects, we have investigated the 
conformations of the network strands as a function of their 
contour length N s . For high q all structure functions fall on 
top of each other and show the same fractal structure with 
S(q) ~ g _1 /°- 72 ( 2 ) (Fig.||). The complementary Fig. |3] shows 
a log-log plot of the mean square strand extension (r 2 ) eq (N s ) 
versus their length. In agreement with the results for the struc- 
ture functions, we find a power law (r 2 ) eq oc fr 2 _/v 2x0 - 72 for 
strands which are shorter than the effective strand length N e 
and therefore sub-affine deformations. Long strands, on the 



other hand, deform affinely with (r 2 ) eq = (r 
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Clearly, the results of our simulations do not agree with the 
predictions of any of the theories presented in the introduc- 
tion. While the neglect of entanglements seems to be fairly 
simple to repair by treating them as effective cross-links (with 
N e supplanting the average strand length N s \ 19]), the fractal 



structure of the strands and the exponent v = 0.72(2) come as 
a surprise. In the following, we discuss a possible explanation 
for the stronger swelling of network strands {y w 0.72) than 
of single chains (v « 3/5) under good solvent conditions. 

We begin by recalling Flory's argument [[!]] for the typical 
size Rp oc bN v of a single polymer chain of length TV and 
statistical segment size b in a good solvent. The equilibrium 
between an elastic energy oc R F /(b 2 N) of a Gaussian chain 
stretched to Rp and a repulsive energy oc b d R F ( N / R F ) 2 due 
to binary contacts between monomers in d dimensions leads 
to 



The simplest models for swollen networks have the regular 
connectivity of a crystal lattice. In agreement with the c*- 
theorem they adopt equilibrium conformations with strand ex- 
tensions of the order of Rp [|28|]. However, these systems are 
hardly good models for the swelling process of networks pre- 
pared in the dry state, since the hypothetical initial state at 
melt density has an unphysical local structure with average 
strand extensions RpQ^ 1 / 3 oc bN^ 3 as in dense globules. 
In contrast, if the corresponding semi-dilute solution is com- 
pressed, the chains shrink only weakly from Rp to the Gaus- 



sian coil radius R 



bN. 
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interpenetrating with n oc N s of them sharing a volume of 
R 3 . Moreover, at least the simplest model for highly cross- 
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linked networks prepared in the dry state, n ~ N s mutu- 
ally interpenetrating regular networks with strand extensions 
of the order of R [24|, cannot possibly comply with the c*- 
theorem, if one disregards macroscopic chain separation: Ei- 
ther the strands extend to Rp, leading to internal concentra- 
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tions of c*N s , or the systems swell to c , in which case 
the strands are stretched to RpN^ 6 . The same conclusions 
should hold for any network without too many defects, where 
the global connectivity forces neighboring chains to share the 
same volume independent of the degree of swelling. 

The consequences can be estimated using a simple Flory 
argument. Instead of a single polymer, we now consider a 
group of chains which can swell but not desinterpenetrate, i.e. 
n ~ N 1 / 2 chains of length N which span a volume R FR . The 
equilibrium between the elastic energy oc n R 2 FR / (b 2 N) and 
the repulsive energy oc b d R FFl (n N/R FR ) 2 leads to 



Instead, they become highly 



d + 2 



(1) 



4 + 2d 



(2) 



Quite interestingly, this local argument reproduces in three 
dimensions with Q eq ~ N d ^ d+2 ^ = N 3 / 5 and R FR ~ 
Q 1 / d bN 1 / 2 ~ TV 7 / 10 the results of the classical Flory-Rehner 
theory of gels. However, in analogy to the Flory argument for 
single chains, Eq. (||) should also apply to subchains of length 
G with 1 G < N which share their volume with a corre- 
spondingly smaller number of other subchains. In particular, 
the local degree of swelling, G 1 / 5 , should be sub-qffine and 
the exponent v — 7/10 should characterize the entire local 
chain structure up to the length scale of the effective strand 
length, N e , This is in excellent agreement with the main find- 
ings from our simulations (see Figs. |2| and ^|). 

Before we conclude, some additional remarks are in order: 
(i) For swelling in a Theta-solvent, the analogous scaling ar- 
gument yields Q ~ TV 3 / 8 in agreement with previous theories 
and experiments [ |l6|]l7| , |T9| ] and predicts local chain structures 
characterized by v = 5/8. (ii) Eq. (||) can also be derived 
along the lines of JT^ , [l9| , p4| ] from an equilibrium between the 
elastic energy of blob chains and the osmotic pressure of a 
semi-dilute polymer solutions. Note, that the apropriate blob 
size is a function of the size G of the subchains under consid- 
eration and that isolated chain behavior is only expected be- 
low the original correlation length ^ prep for systems prepared 
by cross-linking semi-dilute solutions, (iii) Sommer, Vilgis 
and Heinrich [ |29| ] have argued that the effective inner frac- 
tal dimension di of a polymer network is larger than di — 1 
for linear chains, leading to stronger swelling with v = 
While the correction goes into the right direction, it is difficult 
to explain a strand length independent effective inner fractal 
dimension of di = 1.5 as an effect of the local connectivity, 
(iv) However, such effects may well be important in systems 
with a sufficient number of defects such as dangling ends or 
clusters. If the global connectivity is weak, the chains may 
locally desinterpenetrate, leading to a behavior which agrees 
much better with the c* -theorem [M,M]. 

In summary, we have used large scale computer simulations 
and scaling arguments to investigate the swelling behavior of 
defect free model networks prepared at melt density. We find 
that for networks with average strand lengths N s > N e the 
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swelling is limited by entanglements to Q eq ss N e and that 
the strands locally exhibt a fractal structure characterized by 
an exponent v ks 7/10 which should be directly observable in 
neutron scattering experiments. 
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